D  A089668 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison.  Wisconsin  53706 


June  1980 


(Received  June  2,  1980) 


DTIC 

^ELECTE 


Approve#  for  public  roleose 

^  Distribution  unlimited 

£ 

c 

.  Sponsored  by 

uJ 

-1 m  s.  Army  Research  Office 
p .  o.  Box  12211 
Research  Triangle  Park 
^•^Jorth  Carolina  27709 


National  Science  Foundation 
Washington,  D.  C.  20550 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

*  y 

A  NONLINEAR  HYPERBOLIC  VOLTERRA  EQUATION  IN  VISCOELASTICITY*  • 


*  ** 

C.  M.  Dafermos  and  J.  A.  Nohel 


Technical  Summary  Report  #2095 
June  1980 

ABSTRACT 
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A  general  model  for  the  nonlinear  motion  of  a  one  dimensional,  finite, 
homogeneous,  viscoelastic  body  is  developed  and  analysed  by  an  energy 
method.  It  is  shown  that  under  physically  reasonable  conditions  the  nonlinear 
boundary,  initial  value  problem  has  a  unique,  smooth  solution  (global  in 
time),  provided  the  given  data  are  sufficiently  "small"  and  smooth;  moreover, 
the  solution  and  its  derivatives  of  first  and  second  order  decay  to  zero  as 
t  “.  Various  modifications  and  generalizations,  including  two  and  three 
dimensional  problems,  are  also  discussed. 
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SIGNIFICANCE  AND  EXPLANATION 


In  nonlinear  systems  of  "hyperbolic"  type,  characteristic  speeds  are  not 
constant  so  that  weak  waves  are  amplified  and  smooth  solutions  may  blow  up  in 
finite  time  due  to  the  formation  of  shock  waves.  It  is  interesting  to 
consider  situations  where  this  destabilizing  mechanism  coexists  (and  thus 
competes)  with  dissipation. 

In  certain  cases  (e.g.,  viscosity  of  the  rate  type)  dissipation  is  so 
powerful  that  waves  cannot  break  and  solutions  remain  globally  smooth.  A  more 
interesting  situation  arises  when  the  amplification  and  decay  mechanisms  have 
comparable  power  so  that  the  outcome  of  their  confrontation  cannot  be 
predicted  at  the  outset.  Elementary  dimensional  considerations  indicate  that 
breaking  of  waves  develops  on  a  time  scale  inversely  proportional  to  wave 
amplitude  while  dissipation  proceeds  at  a  roughly  constant  time  scale.  It 
should  thus  be  expected  that  dissipation  prevails  and  waves  do  not  break  when 
the  initial  data  are  "small".  Results  of  this  type  were  obtained  by  T.  Nishida 
for  the  guasilinear  wave  equation  with  first-order  frictional  damping  for 
sufficiently  smooth  and  small  initial  displacements  and  initial  velocities. 

In  this  paper  we  develop  and  study  a  general  nonlinear  model  for  the 
motion  of  a  one  dimensional,  finite,  homogeneous  body.  Here  the  dissipation 
mechanism  which  is  induced  by  memory  effects  of  the  viscoelastic  materials 
( stress-strain  relaxation  function  -  the  stress  is  a  nonlinear  functional 
rather  than  a  function  of  the  strain)  is  different  and  more  subtle.  Using 
elementary  energy  methods,  which  are  combined  with  frequency  domain  techniques 
for  nonlinear  Volterra  equations,  we  show  that  under  physically  reasonable 
conditions  on  the  stress-strain  relaxation  function,  the  known  history  of  the 
displacement,  the  nonlinearities  of  the  model,  and  on  the  assigned  external 
body  force,  the  boundary-history  value  problem  (1.9),  (1.18)  in  the  text  which 
describes  the  model  has  a  unique,  smooth  solution  (global  in  time),  provided 
the  given  data  (history  and  external  body  force)  are  sufficiently  smooth  and 
"small".  Moreover,  we  also  show  that  the  solution  and  its  spacial  and  time 
derivatives  of  first  and  second  order  decay  to  zero  as  t  “.  Various 
modifications  and  generalizations  of  the  model,  including  two  and  three 
dimensional  problems,  are  also  considered. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  NONLINEAR  HYPERBOLIC  VOLTERRA  EQUATION  IN  VISCOELASTICITY 

*  ** 

C.  M.  Dafermos  and  J.  A.  Nohel 

1.  Introduction.  In  nonlinear  systems  of  "hyperbolic"  type,  characteristic  speeds 
are  not  constant  so  that  weak  waves  are  amplified  and  smooth  solutions  may  blow  up 
in  finite  time  due  to  the  formation  of  shock  waves.  It  would  be  interesting  to 
consider  situations  where  this  destabilizing  mechanism  coexists  (and  thus  competes) 
with  dissipation. 

In  certain  cases  (e.g.,  viscosity  of  the  rate  type)  dissipation  is  so  powerful 
that  waves  cannot  break  and  solutions  remain  globally  smooth.  A  more  interesting 
situation  arises  when  the  amplification  and  decay  mechanisms  have  comparable  power 
so  that  the  outcome  of  their  confrontation  cannot  be  predicted  at  the  outset. 
Elementary  dimensional  considerations  indicate  that  breaking  of  waves  develops  on  a 
time  scale  inversely  proportional  to  wave  amplitude  while  dissipation  proceeds  at  a 
roughly  constant  time  scale.  It  should  thus  be  expected  that  dissipation  prevails 
and  waves  do  not  break  when  the  initial  data  are  "small".  Results  of  this  type  for 
quasilinear  wave  equations  with  frictional  damping  were  first  obtained  by  Nishida 
[1]  and,  subsequently,  by  Matsumura  [2],  who  uses  methodology  that  goes  back  to 
Schauder  (3].  The  more  delicate  situation  of  thermal  damping  (one  dimensional 
thermoelasticity)  is  discussed  in  Slemrod  [4]. 

A  different,  subtler  type  of  dissipation  mechanism  is  induced  by  memory  effects 
and  arises  in  nonlinear  viscoelasticity.  A  simple,  one  dimensional,  model 
corresponds  to  the  constitutive  relation 

t 

(1.1)  o(t,x)  *  v?(e(t,x))  +  /  a'(t  -  T)i|>(e(T,x))di  , 
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where  0  is  the  stress,  e  the  strain,  a  the  relaxation  function  with  '  =  d/dt, 
and  v?,  ^  assigned  constitutive  functions.  We  normalize  the  relaxation  function  so 
that  a(")  =  0.  When  the  reference  configuration  is  a  natural  state, 
v'(0)  =  i^(0)  =  0.  Experience  indicates  that  s?(e),  <Me),  as  well  as  the 
equilibrium  stress 

dc  f 

(1.2)  X(e)  =  <^(e)  -  a(0)ty(e) 

are  increasing  functions  of  e,  at  least  near  equilibrium  (e  small).  Moreover,  the 

effect  of  viscosity  is  dissipative.  To  express  mathematically  the  above  physical 

requirements,  we  impose  upon  a(t),  <£(e),  \J>(e)  and  \(e)  the  following  assumptions: 
2  1 

(1.3)  a(t)  e  W  '  (O,00),  a(t)  is  strongly  positive  definite  on  (0,“>)  ; 

(1.4)  v?(e)  e  C3(-“,«),  ^(0)  =  0,  <e’(0)  >  0  ; 

(1.5)  <)(e)  6  C3(-«,<»),  iHO)  =  0,  ^'(0)  >  0  ; 

(1.6)  X’(0)  =  ^'(0)  -  a(0)i|>'(0)  >  0  . 

Assumption  (1.3),  which  requires  that  a(t)  -  a  exp(-t)  be  a  positive  definite 
kernel  on  [0,«)  for  some  a  >  0,  expresses  the  dissipative  character  of 
viscosity.  Smooth,  integrable,  nonincreasing,  convex  relaxation  functions,  e.g., 

K 

(1.7)  a(t)  =  £  \>kexp<-ukt),  \  >  °'  \  >  0  * 

which  are  commonly  employed  in  the  applications  of  the  theory  of  viscoelasticity, 
satisfy  (1.3). 

It  is  often  convenient  to  express  a,  given  by  (1.1),  in  terms  of  equilibrium 
stress,  namely  (integrate  (1.1)  by  parts  and  use  (1.2)), 

t 

(1.8)  o(t,x)  =  x  ( e(  t  ,x) )  +  /  a(t  -  T)ip(e(T,x))  dx  . 

—CO  ^ 

We  now  consider  a  homogeneous,  one  dimensional  body  (string  or  bar)  with 
reference  configuration  [0,1]  of  density  p  =  1  (for  simplicity)  and  constitutive 
relation  (1.1)  which  is  moving  under  the  action  of  an  assigned  body  force  g(t,x), 
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-00  <  t  <  « 


04x41.  We  let  u(t,x)  denote  the  displacement  of  particle  x  at 


time  t  in  which  case  the  strain  is  e(t,x)  ■  ux(t,x).  Thus  the  equation  of  motion 

pu  ^  =o  +  pg  here  takes  the  form 
tt  x  ’ 

t 

(1.9)  u  =  >p(u  )  +  /  a'(t  -  T)ij/(u  )  dr  +  g,  -“  <  t  <  “>,  0  4x41, 

ttxx'  XX 

«Q0 

or,  if  one  uses  representation  (1.8)  for  the  stress, 

t 

(1.10)  u  =  x(u  )  +  /  a(t  -  T)iJ>(u  )  dt  +  g,  -“  <  t  <  ~,  0  4x41. 

tt  X  X  X  TX 

■OP 

The  history  of  the  motion  of  the  body  up  to  time  t  =  0  is  assumed  known,  i.e., 

(1.11)  u(t,x)  ■  v(t,x),  -*  <  t  4  0,  0  4x41  , 

where  v(t,x)  is  a  given  function  which  satisfies  equation  (1.9)  together  with 
appropriate  boundary  conditions,  for  t  <  0.  Our  task  is  to  determine  a  smooth 
extension  u(t,x)  of  v(t,x)  on  (-*,*)  x  [0,1]  which  satisfies  (1.9)  together 
with  assigned  boundary  conditions,  for  -«•  <  t  <  ». 

Upon  setting 

0 

(1.12)  h  =  J  a'(t  -  x )i|/(v  )  dt  +  g,  t  >  0,  0  4x41, 

M  XX 

■w 

(1.13)  uQ(x)  *  v(0,x),  u^x)  =  vt(0,x),  0  4x41  , 

the  history-value  problem  (1.9),  (1.11)  reduces  to  the  initial-value  problem 

t 

(1.14)  u  *v?(u  )  +  /  a'  ( t  -  T )  \|i(  u  )  dT  +  h,  0  4  t  <  »,  0  4x41, 

tt  x  x  xx 

(1.15)  u(0,x)  -  uQ(x),  ufc(0,x)  -  u,(x),  0<x<1  . 

Conversely,  (1.14),  (1.15)  can  be  reduced  to  (1.9),  (1.11)  by  constructing  a  func¬ 
tion  v(t,x)  on  (—,0)  x  (0,1)  which  satisfies  v(0,x)  ■  uQ(x),  vt(0,x)  =  u^x), 


3 


vtt(0,x)-*<u0x)x+h(0,x)  . 


(1.16)< 


v  (0,x)«^"(u  )u  u  +<fi' ( u  )u  +a'(0)<i(uA  )  +h  (0,x), 
ttt'  Ox  Oxx  lx  Ox  Ixx  Ox  x  t 


0  <  x  <  1, 
0  <  x  <  1, 


together  with  appropriate  boundary  conditions,  for  t  <  0,  and  then  defining 
g(t,x)  on  (-“,«“)  x  [0,1]  by 


v. 

f  V..  -  I p{v  )  -  /  a'(t-t)ij/(v  )  dT, 

I  tt  V  V  J  ~  ~ 


(1.17)  g(t,x)  -< 


X  X 


X  X 


h  -  /  a' (t-T )i|)(\’x)xdT , 


t  <  0,  0  <  x  <  1  , 

t  >  0,  0  <  x  <  1  . 


The  purpose  of  (1.16)  is  to  ensure  that  g(t,x),  as  defined  by  (1.17),  has  the 
smoothness  properties,  across  t  *  0,  which  will  be  required  below  in  the  existence 
theorem. 

For  the  special  case  i>(e)  =  <fi(e)  variants  of  existence  theorems  for  (1.14), 
(1.15)  were  established  by  MacCamy  [5],  Dafermos  and  Nohel  [6]  and  Staffans  [7], 

The  assumption  i|i  =  allows  one  to  invert  the  linear  Volterra  integral  operator  on 
the  right-hand  side  of  (1.14)  and  thus  express  ^ ( ux ) x  in  terms  of  utt  -  h 
through  an  inverse  Volterra  integral  operator  using  the  resolvent  kernel  associated 
with  a*.  One  may  then  transfer  time  derivatives  from  ufct  to  the  resolvent  kernel 
via  integration  by  parts.  This  procedure  reveals  the  instantaneous  character  of 
dissipation  and,  at  the  same  time,  renders  the  memory  term  linear  and  milder,  thus 
simplifying  the  analysis  considerably.  On  the  other  hand,  the  above  approach  is 
somewhat  artificial:  By  inverting  the  right-hand  side  of  (1.14),  one  loses  sight  of 
the  original  equation  and  of  the  physical  Interpretation  of  the  derived  a  priori 
estimates.  More  importantly,  the  physical  appropriateness  of  the  restriction 
ii  *  is  by  no  means  clear. 

Remark:  The  present  normalisation  of  the  kernel  a  with  a(®)  *  0  is  different 
from  that  in  the  existing  literature  (see  [5],  [6],  [7]).  The  reader  should  note 
a',  not  a,  enters  the  constitutive  relation  (1.1)  as  well  as  the  equation  of 


•4 


motion  (1.9).  The  present  normalization  is  more  convenient  for  technical  reasons; 

for,  the  equivalent  form  (1.8)  of  the  constitutive  equation  in  which  a  (rather 

than  a')  enters,  and  the  corresponding  equation  of  motion  (1.10),  are  extremely 

convenient  for  obtaining  the  crucial  a  priori  estimates  in  our  analysis,  when  a 

satisfies  assumption  (1.3).  In  the  earlier  literature  in  which  only  the  special 

case  =  &  was  studied,  the  normalization 

a(t)  ■  a  +  A(t)  0  <  t  <  •  , 

2  1 

a(0)  ■1,a##>0,  AeW*  (0,*),  A  strongly  positive  was  used. 

In  another  noteworthy  special  case,  when  aft)  *  exp(-yt),  (1.9)  is  equivalent 
to  the  third  order  partial  differential  equation 

ut«  *  +  "’"V*  *  «t  ♦  «« 

studied  by  Greenberg  [8] . 

In  this  paper  we  show  how  one  may  deal  with  Equation  (1.9)  directly  and 
establish  existence  of  solutions  without  the  assumption  V  =  i|>.  We  will  consider  in 
detail  the  case  where  the  boundary  of  the  body  is  free  of  traction  which  leads  to 
boundary  conditions 

(1.18)  0(t,O)  ■  0(t,1)  -  0,  -«  <  t  <  •  . 

Other  types  of  boundary  conditions  will  be  discussed  in  Section  4.  The  change  of 
variable  (superposition  of  a  rigid  motion) 


t  t  1 


(1.19) 


i(t,x)  *  u(t,x)  +  m  +  m  t  +  /  /  /  g(s,y)dydsdr 
0  1  0  0  0 


shows  that  without  loss  of  generality  we  may  assume 


(1.20) 


1 

/  g(t,x)dx  *  0, 
0 


<  t  <  •  , 


(1.21) 


/  u(t,x)dx  *  0, 
0 


<  t  < 
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Of  the  body  force  g  we  require 


(1.22) 


|g(t,»),gt(t,»),gx(t,«)  in  C((-“,a>)f  L2(0,1))n  L2((-“,“>);  L2(0,1)) 
[g(t,x)=g,j(t,x)+g  (t,x)  with  g1tfc(  t ,  • )  »g2tx< t '  *  >  in  l2^  <  jL2(0 , 1 ) ) . 


As  noted  above,  despite  the  presence  of  viscous  dissipation,  it  is  not  to  be 
expected  that  a  global  smooth  solution  to  (1.9),  (1.11),  (1.18)  will  exist  unless 
the  amplitude  of  waves  remains  small.  Consequently,  one  may  only  hope  to  obtain 
global  existence  results  under  the  restriction  that  g(t,x)  be  appropriately 
"small".  We  "measure"  g(t,x)  by 


(1.23)  G  d=f  sup  /  {g^+g^+gx)(t,x)dx  +  /  /  {g^+g^+g^+g^+g^Jdxdt 

(-00,00)  0  0  U 


1 


2.  2.  2- 


00  ) 


2.  2.  2.  2 


Our  main  result  is 

Theorem  1.1.  Under  assumptions  (1.3),  (1.4),  (1.5),  (1.6),  there  exists  a  constant 
V  >  0  with  the  following  property:  For  every  g(t,x)  £n  (-“,“)  x  [0,1]  which 
satisfies  (1.20)  and  (1.22)  with 

(1.24)  G  <  v2  , 

and  for  any  v(t,x)  on  (-<",0]  x  [0,1],  with  v(t,*),  v^ft,*),  v^(t ,•),  v^ft,*), 

vtx(t'#)'  vxx(t'°'  Vttt(t'*)'  vttx(t'*5'  Vtxx(t'°'  Vxxx(t'0  -is. 

C((-°°,0];  L2(0,D)  O  L2((-°°,0];  L2(0,1)),  which  satisfies  Equation  (1.9)  together 

with  the  boundary  conditions  (1.18)  for  t  <  0,  there  exists  a  unique  u(t,x)  on 

(_oo,oo)  x  [0,1],  with  u  ( t ,  • ) ,  u  (t,*),  u  (t,»),  u  (t,»),  u  (t,«),  u  ( t ,  • ) , 

u  X  lu  lX  XX 

Uttt(t',}'  Uttx(t'#)'  Utxx(t'*)'  Uxxx(t'#)  is. 

C( ( -*,*) ;  L2(0,1))  n  L2 ((-“,»);  L2(0,1)),  which  satisfies  (1.9),  (1.11),  (1.18), 
as  well  as  (1.21).  Furthermore, 

(1.25)  u(t,*),ut(t,*),ux(t,»),utt(t,»),utx(t,«),uxx(t,*)  SSllj*  0,  t  ®. 

The  proof  of  the  above  theorem  employs  the  general  strategy  developed  in 

[2,6,7].  We  first  establish,  in  Section  2,  the  existence  of  a  local  solution. 


defined  on  a  maximal  interval  (-“,TQ),  with  the  property  that  when  TQ  <  ®  a 
certain  norm  blows  up  as  t  +  T  .  Then,  in  Section  3,  we  show  that,  due  to  viscous 
dissipation,  the  afforementioned  norm  remains  uniformly  bounded  on  the  maximal 
interval,  provided  that  (1.24)  holds  with  w  sufficiently  small.  In  particular, 

TQ  =  »  and  the  smooth  solution  exists  globally. 

In  the  final  Section  4,  we  have  collected  information  on  various  extensions  of 
the  above  results.  We  show  how  one  can  handle  boundary  conditions  other  than 
(1.18).  We  indicate  how  alternative  sets  of  assumptions  on  v(t,x)  and  g(t,x) 
lead  to  variants  of  Theorem  1.1  rendering  information  on  the  smoothness  of 
solutions.  Finally,  we  explain  how  the  present  techniques  may  be  used  to  establish 
existence  theorems  for  the  equations  of  multidimensional  viscoelasticity  as  well  as 
abstract  integrodif ferential  equation  in  Hilbert  space. 

2.  Local  Solutions.  In  this  section  we  establish  a  local  existence  theorem  on  a 
maximal  interval.  It  is  more  convenient  to  work  with  Equation  (1.14)  to  which,  as 
we  have  seen,  (1.9)  may  be  reduced.  Also  we  shall  impose  here  boundary  conditions 

(2.1)  ux(t,0)  =  u^t.l)  -  0,  t  >  0  , 

which,  though  apparently  stronger  than  (1.18),  are  actually  equivalent  to  (1.18),  as 
will  be  shown  in  Section  3.  Finally,  we  temporarily  strengthen  assumption  (1.4) 
into 

(2.2)  i£(  e)  e  C3(-“,»),  ^(0)  =  0,  *’(e)  >  x  >  0,  -“>  <  e  <  “>  . 

On  the  other  hand,  assumptions  ^’(O)  >  0,  >  0  and  the  positivity  of  the 

kernel  a(t)  will  not  play  any  role  in  this  section. 

Theorem  2.1.  Let  u^  (x),  u0jf(x),  u0xx(x),  Uoxxx(x) '  u1(x>'  u1x(x)' 
uixx<x>  be  in  L2(0,1)  and  assume 

(2.3)  u0x(0)  =  Uqx(  1 )  =  0,  u1x(0)  *  uixM)  *  0  • 

Moreover,  let  h(t,x)  be  defined  on  [0  ,<*»)  x  (0,1]  with  h(t,»),  h^(t,*),  h_(t,*) 
in  C([0,«);  L2(0,1))  and  h(t,x)  *  ht(t,x)  +  h2(t,x),  h1tt,  h2tx  in. 
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2  2 

L  (10,-);  L  (0,1)).  Then  there  is  TQ,  0  <  T^  <  ®,  and  a  unique  function 
u ( t , x)  e  C2((0,Tq)  x  [0,1)),  with  u  uttx(t,*),  Uxxx<t#*5  — 

C( [0,T] j  L2(0,1)),  for  every  0  <  T  <  TQ,  such  that  u  satisfies  (1.14)  on 
[0 ,Tg )  x  [0,1]  together  with  initial  conditions  (1.15)  and  boundary  conditions 
(2.1)  on  [0,Tq).  Furthermore,  if  TQ  <  », 

1  2  o  2  2  2  2  2 

/  {u  (t,x)+u  (t,x)+u  (t,x)+u  (t,x)+u  (t,x)+u  (t,x)+u  (t ,x) 

(2.4)  0 

+  Uttx(t'x)+Utxx(t'x)+Uxxx(t'x)}dx  *  +  T0  • 

We  note  that  h(t,x)  and  Ug(x),  u^(x),  defined  by  (1.12),  (1.13)  with 
v(t,x)  and  g(t,x)  as  in  Theorem  1.1,  do  satisfy  the  assumptions  of  Theorem  2.1. 

The  proof  of  Theorem  2.1  which  is  a  variant  of  the  local  result  in  [6]  will  be 
based  upon  the  Banach  fixed  point  theorem.  We  begin  with  some  preparation.  For 
M,T  >0,  we  let  X(M,T)  denote  the  set  of  functions  w(t,x)  on  [0,T]  x  [0,1], 

with  w( t , • ) ,  w  (t,«),  wx(t,*),  wfct(t,«>,  wxx<t'*>'  Wtttft'*>'  Wttxft'*>' 

w  (t,*),  w  (t,*)  in  L*([0,T);  L2(0,1))  which  assume  initial  data 
txx  XXX 

w(0,x)  =  uQ(x),  wfc(0,x)  =  u,j(x)  and  boundary  conditions  wx(t,0)  =  wx(t,1)  =  0, 
t  e  [0 ,T] ,  and  satisfy 

_  ^  A  .  .  .  . 

(2.5)  ess-sup  /  {w  (t,x)  +  w  (t,x)  +  w  (t,x)  +  w  (t,x)}dx  <  M  . 

[0  ,T]  0  ttX  tXX  XXX 

For  w(t,x)  e  X(M,T) ,  (2.5)  and  the  Poincare  inequality  yield 

(2.6)  w2(t,x)  +  w2x(t,x)  +  w^(t,x)  <  M2,  0  <  t  <  T,  0  <  x  <  1  . 

2 

We  now  consider  the  map  S  :  X(M,T)  +  C  ([0,T]  x  [0,1])  which  carries 
w(t,x)  e  X (M,T)  into  the  solution  u(t,x)  of  the  linear  equation 

t 

(2.7)  u  -^’(w  ) u  =  /  a'(t  -  T)<Mw  )  dx  +  h 

tt  x  xx  o  xx 

satisfying  initial  conditions  (1.15)  and  boundary  conditions  (2.1).  we  note 


■8- 


are  in 


that  <fi ' <wx(t,x) )  is  C1  smooth  and  v5'(wx>tt,  ^’(wx>tx 
00  2 

L  ( (0 ,T] ;  L  (0,1)).  Furthermore,  if  f(t,x)  denotes  the  right-hand  side  of  (2.7) 
then  f(t,*),  ft(t,«),  fx(t,»)  are  in  C(tO,T)>  L2(0,1))  and  f(t,x)  =  f.,<t,x) 

f2(t,x)  with  f1tt,  f2tx  in  ^2(tO,T]>  L2  (0 , 1 ) ) •  It  then  follows  by  standard 

theory  that  uttt<t,«),  uttx(t,,)'  utxx(t,,)  and  uxxx(t,,)  are  in  c( f°'T5 ; 
L£(0,1)).  Our  strategy  is  to  show  that,  under  proper  conditions,  S  has  a  unique 
fixed  point  in  X(M,T)  which  will  obviously  be  the  solution  to  (1.14),  (1.15), 
(2.1)  with  the  desired  properties. 

Lemma  2.1.  When  M  is  sufficiently  large  and  T  is  sufficiently  small,  S  maps 
X(M,T)  into  itself. 

Proof.  We  fix  n  >  0  and  apply  to  (2.7)  the  forward  difference  operator  A, 
(Aw)(t)  «  u(t  +  n)  -  u>(t) ,  thus  obtaining 


(2.8) 


Au  .  -  *  (w  )Au 

tt  X  XX 


A^  *  (w  )u  +  Aip*  (w  )Au 

Y  X  XX  X  XX 

t 

+  A  /  a'(t-T)<J>(w)dT+Ah. 

0  X  x 


We  multiply  (2.8)  by  ^utxx  and  integrate  over  [0,s]  x  [0,1],  0  <  s  <  T.  After 

2 

appropriate  integrations  by  parts,  we  divide  through  by  n  and  we  let  n  +  0.  We 
give  the  details  of  the  computation  of  one  term: 


(2.9) 


1 


I  f 
0  0 


Au  Au  dxdt 
tt  txx 


s  1 

/  I 
0  0 


u 


Au^  Au  dxdt 
ttx  tx 


( Au 


tx 


y  1 

)  (s,x)dx  +  y  / 


(Aii^)  ( 0 ,  x )  dx 


whence 


1 

(2.10)  lim  —  /  /  Au  Aufcxxdxdt 
n+0  n  0  0 


where 


-9- 


(2.11) 


utt x<0'x)  *  *(uOx(x>)xx  +  V°'x)  * 


We  apply  the  same  procedure  to  the  remaining  terms  of  (2.8)  thus  obtaining 


(2.12) 


1  i  -  .  1  ,  -  1 

2  /  uHJs'xl(fe  +  7/  *'<w  (s,x))u  (s,x)dx  =  X  /  u  (0,x)dx 


2  '  “ttx 


2  '  “ttx' 


1  1  ?  1  S  1  2 
+  2  /  (U0x(X>,U1XX(X)dX  +  2  /  / 


x  tx  txx 


si  si 

+  /  /  * (W  )w  w  u  u  dxdt  +//  ifi"(w  )w  u  u  dxdt 

q  q  x  tx  xx  xx  ttx  q  q  x  txx  xx  ttx 

s  1  1 

+  /  /  >f"(w  )w  u  u  dxdt  -  /  a'(s)iHu„  (x))  u  (s,x)dx 

q  q  x  tx  xxx  ttx  '  Ox  x  txx 


1  s  1 

+  /  a'(0)ij>(u  (x))  u.  (x)dx  +  /  /  a"(t)f(u„  )  u  dxdt 

0  Ox  x  Ixx  '  *  Ox  x  txx 


/  u*xJs'x)  /  a'(s  -  t)\|i(w  <t,x))  dtdx  +  J  /  a'(0)t(w  )  u  dxdt 
0  0  xx  0  0  XX  txx 


+  /  /  u  /  a"(t  -  T)\p(w  )  dTdxdt  +  /  /  h  u  dxdt 
f  4  fvv  4  v  v  44  2tX  ttx 


J  J  txx  ' 

oo  tx  o 


1  1  s  1 

-  /  h  (s,x)u  (s,x)dx  +  /  h  (0,x)u  (x)dx  +  /  /  h„  u  dxdt  . 
q  It  txx  '  It  Ixx  q  *  Itt  txx 


We  now  differentiate  (2.7)  with  respect  to  t  and  x  to  obtain 


(2.13) 


u  -^'(w  )u  -  ^"(u.  )u„  u,  -  /  fa"(W  )w  u  }  dT 
ttt  X  txx  Ox  Oxx  lx  J0  X  TX  XX  T 


a'(t)i^(u  )  +  f  a'(t  -  t  )  ^  { w  )  dt  +  h  , 

Oxx  '  X  TX  t 
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(2.14) 


U  ^  -  *?'  (w  )U  -  S?"(u„  )U-  -  /  {^"(w  )w  U  }  dt 

ttx  X  xxx  Ox  Oxx  J  X  XX  XX  T 


a*  (t  -  T  )tjf(w  )  dT  +  h 
X  XX  X 


from  which  we  easily  get  the  estimates 


(2.15) 


1  2  1  2  2 
/  u  ..  (s,x)dx  -  6  /  (s,x))  u.  (s,x)dx 

t  ttt  i  x  txx 

0  0 


1 


<6/  <fi " ( uQx( x ) ) 2Upxx( x ) u^x( x ) dx 


1  s  1  ,  ? 

+  6  /  {/  [*"(w  >w  u  ]  dt}  dx  +  6  /  a'(t)  i|;(u  (x))  dx 

00  x  tx  xx  r  Q  ux  x 


Is  1  s  2 

+  6  /  {/  a*(s  -  t)ij>(w  )  dt)  dx  +  6  /  /  h  dxdt  , 

0  0  x  tx  oo 


(2.16)  /  ^(wx(s>x))2u^x(s,x)dx  -  5  /  u2tx(s,x)dx 


1  .2  4  18 


<  5  /  *"(u0x(x>)  u0xxdx  +  5  /  {/  t'r(wx)WxxUxx1tdt}  dX 


0  0 


Is  1  s 

+  5  /  {/  a'(s  -  t)ii»(w  )  dt}  dx  +  5  /  /  h  dxdt  . 

0  0  *  xx  0  0  x 


Let  us  set 


1 


2  f '  f  2  2  2  2  ^  2  ^  2  ^  2 

N  ‘l  {u0  +  U0x  *  U0xx  +  U0xxx  +  U1  +  U1x  +  U1xx}dX 


(2.17) 


1  <*>  1 
+  sup  /  {h2(t,x)+h2(t,x)+h2(t,x)}dx  +  /  /  ^tt+h2tx^dxdt  ' 

(0,“)  0  x  0  o 
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Then,  by  virtue  of  (2.5),  (2.6),  the  Poincare  inequality  and  Schwarz's  inequality, 
every  term  on  the  right-hand  side  of  (2.12),  (2.15),  and  (2.16)  can  be  majorized  by 
one  of  p(N) ,  Tq(M)V2,  p(N)V,  Tq(M)V,  T2q(M)V2,  T^p(N)V,  T2q(M),  Tp(N),  where 
p(»)  and  q( • )  are  locally  bounded  functions  on  [0,°®).  Thus,  combining  (2.12), 
(2.15)  and  (2.16)  and  using  (2.2),  we  arrive  at  an  estimate  of  the  form 

(2.19)  V2  <  c{p(N)  +  Tq(M)V2  +  p(N)V  +  Tq(M)V  +  T2q(M)V2 

+  T1/2p(N)V  +  T2q(M)  +  Tp(N)}  . 

Applying  the  Cauchy-Schwarz  inequality, 

(2.20)  {1  -  cTq(M)  -  2cT2q(M)  -  cTp(N)  -  ^ }V2 

<  c{2p(N)  +  cp2 (N)  +  q(M)  +  T2q(M)  +  Tp(N)}  . 

2  2 

Thus,  if  one  fixes  M  >  8c{2p(N)  +  cp  (N)  +  q(M)}  and  then  selects  T  so  small 

2 

that,  at  the  same  time,  cTq(M)  +  2cT  q(M)  +  CTp(N)  <  %  and 
cT2q(M)  +  cTp(N)  <  M2/8,  (2.20)  yields  V2  <  M2  and  u(t,x)  e  X(M,T).  The  proof 
of  the  lemma  is  complete. 

We  now  equip  X(M,T)  with  the  metric 


(2.21)  p(u,u) 


(u  -u  )2] (t,x)dx}  ^ 

XX  XX 


where  u,  u  e  X(M,T).  On  account  of  the  lower  semicontinuity  property  of  norms  in 
Banach  space,  X(M,T)  is  complete  under  p. 


Lemma  2.2.  For  M  sufficiently  larqe  and  T  sufficiently  small,  the  ma; 


Proof.  Let  w(t,x),  w(t,x)  e  X(M,T).  We  set  u  =  Sw,  u  =  Sw,  W  =  w  -  w. 


U  =  u  -  u.  Then  U(t,x)  is  the  solution  of  the  problem: 


(2.22) 


where 


°tt  ‘  *'<WxIUxx  *  l  ”txaT 


+  /  a'  (t  -  t)[*'(w  )W  +  B<T'x>«xxwxldT  ' 


(2.23) 


U(0,x)  =  0,  ut(0,x)  =  0,  0  <  X  <  1 


(2.24) 


Ux(t,0)  =  Ux<t,D  =  0,  0  <  t  <  T 


(2.25)  A( t  »x)  = 


(t,x))  -^'(w  (t,x)) 

- 2L. - 2 -  ,  wx(t,x)  *  wx(t,x) 

w  (t,x)  -  w  (t,x) 

X  x 


tfi  "(w^tt.x) ) , 


w  (t,x)  =  w  (t,x)  , 

X  X 


(2.26)  B(t,x)  = 


\|»*  (w  (t|X) )  -  'll'  (w  (t,x) ) 
x  x 


w  (t,x)  -  w  (t,x) 
X  x 


w  (t,x)  *  w  (t,x) 

X  * 


\j)"(wx(t,x) )  > 


w  (t,x)  =  W  (t,x)  . 
X  x 


Furthermore, 


(2.27) 


-  *’<VUtxx  =  +  ^”(Wx)WtxUxx  +  ^"(wx)UxxWtx 


C( t,x)w  W  +  a'(0)f  (*  )W  +  a' (0)B(t,x)w  w 


t  .fc  _ 

+  f  a" ( t  -  T)^' (w  )W  dT  +  J  a"(t  -  T)B(T,x)w  W  dT 

0  x  xx  0 


where 
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(2.28)  C(t,x) 


y'"(w  ( t  /  X )  )  -  <pM(w  (t,x)) 

X _ X _ 

w  (t,x)  -  w  (t,x) 

X  x 


w  (t,x)  *  W  (  t  ,  X  ,  , 

X  X 


sS  (w^(t,x) )  , 


w  (t,x)  =  w  (t,x) 

X  X 


Multiplying  Equation  (2.27)  by  U^t  and  integrating  over  [0,1]  x  [0,s], 
0  <  s  <  T,  we  obtain,  after  certain  integrations  by  parts, 


(2.29) 


7  I  (s,x)dx  +  -r  /  *>' <w  (s,x))uj  <s,x)dx 


2  1  tt 
0 


1  s  1  si 

=  T  /  /  V'"(WX)WtXUtXdXdt  *  J  /  ^"(Wv)WvvUttUtvdXdt 


X  XX  tt  tx 


si  si 

+  /  /  Ru  W  U  dxdt  +  /  /  i^"(w  )w  U  U  dxdt 

i  J  txx  x  tt  i  :  x  tx  tt  xx 

0  0  0  0 


si  si 

+  /  /  ^"(W  )u  V*  U  dxdt  -  /  f  cC;  W  U  dxdt 

i  i  x  xx  tx  tt  ^  tx  x  tt 

0  0  0  0 


si  si 

+  a’(0)  /  f  f (w  )W  U  dxdt  +  a'(0)  f  /  Bw  W  U  dxdt 
0  0  x  xx  fct  00  xx  x  tfc 


si  t 

+  f  f  U  f  a"(t  -  x)<I)'(w  )W  dxdxdt 

o  0  fct  0  x  xx 


si  t 

+  [  [  U  /  a"(t  -  x)Bw  w  dxdxdt  . 
0  0^0 


Moreover,  from  (2.22)  we  get 


• 1  2-2  .  rS  ..  .  .2 


/  '  ( w  (s,x))u  ( s ,  x )  dx  <  3  /  U  (s,x)dx+3  /  Au  (/  W  dx)  dx 

x  xx  0  t  0  XX  0 


Is  2 

+  3  /  {/  a'(s  -  t)[^'(w  )W  +  Bw  W  ]dt}  dx  . 

1  i  x  xx  xx  x 

0  0 


Combining  (2.29)  with  (2.30)  and  using  (2.2),  (2.5),  (2.6),  the  Poincare  inequality 


and  the  Cauchy-Schwarz  inequality,  we  arrive,  a iter  a  long  computation,  at  an 
estimate  of  the  form 


(2.31) 


/  {U  (s,x)  +  U  (s,x)  +  U  <s,x)}dx 
q  tt  tx  xx 


<  (T  +  T2)  max  /  {W2  (t,x)  +  W2  (t,x)  +  W2  (t,x)}dx 
(0  ,T]  0  tt:  tX 

s  1 

+  m  J  J  {U  (t,x)  +  U  (t,x)  +  U  (t,x)}dxdt  , 

0  0  fct  tx  xx 


where  m  depends  solely  upon  a(t),  M,  and  bounds  of  <p ,  i|>  and  their  derivatives 


on  the  interval  [-M,M] .  In  order  to  assist  the  reader  to  see  how  (2.31)  is  derived 
from  (2.29),  (2.30),  we  give  the  details  of  the  estimation  of  one  of  the  most 


complicated  terms  on  the  right-hand  side  of  (2.29): 


(2.32) 


si  t 

/  /  U  J  a"  (t  -  t)i|»'(w  )W  dTdxdt 
0  0  tfc  0  x  xx 


s  1 

2  r  r  ..2 


1  t 


<  7  J  J  U  dxdt  +  2e  /  /  {J  a”(t  -  t)\|>’(w  )w  dt}  dxdt 

e  o  o  “  000  xxx 


<  -  /  /  V  dxdt 
0  0 


S  t 

+  2e  /  {/  ]  a "  ( t 

0  0 

<  t  H  Uttdxdt 
0  0 


1  1  ?  ? 
t ) | dT } { /  |a"(t  -  T)|  /  i|i'(w  )  W  df}dxdt 

o  0  x  ** 


+  2es{/  |  a"  (  t  )  |  dr } 2  {  max  <|>’(e)2}{  sup  /  w2  (t,x)dx}  . 

0  [ 0 ,T)  0  XX 


From  (2.31)  and  Gronwall’s  inequality  we  deduce 


\ 
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(2.33) 


1  2  2  p 

max  /  {u.  (t,x)  +  U  (t,x)  +  U  (t,x)}dx 

„  tt  tx  xx 

[0,T]  0 


<  (T  +  T^)emT  max  /  {wf  (t,x)  +  (t,x)  +  (t,x)}dx  . 

[0  ,T]  0  fct  tX 


2  1 

Thus,  when  T  is  so  small  that  (T  +  T  )exp(mT)  <  —  ,  (2.33)  yields 

1 


p ( Sw , Sw )  <  —  p(w,w),  for  w,w  e  X(M,T) 


(2.34) 

and  the  proof  of  the  lemma  is  complete. 

Proof  of  Theorem  2.1.  From  Lemma  2.2  and  the  Banach  fixed  point  theorem  we  deduce 

the  existence  of  a  unique  fixed  point  of  S  in  X(M,T),  for  conveniently  large 

M  and  appropriately  small  T,  which  will  be  the  unique  solution  of  (1.14),  (1.15), 

(2.1)  on  [0,T]  x  [0,1].  Let  T^  <  00  be  the  maximal  interval  of  existence  of  a 

solution  u ( t , x )  to  (1.14),  (1.15),  (2.1)  with  ut(t,»),  u^(t,*),  u  (t,*), 

u  (t,«),  u  (t,*),  u  (t,*),  u  (t,»),  u  ( t ,  • ) ,  u  (t,*)  in 

tx  XX  ttt  ttx  txx  XXX 

00  2 

L  ( [ 0 , T] ;  L  (0,1))  for  every  0  <  T  <  Tfl .  If  Tq  <  “  and  (2.4)  is  not  satisfied, 

2 

we  can  extend  u(t,x)  up  to  t  =  Tg  so  that  u(t,x)  C  ([0,TQ]  x  [0,1]). 

Moreover,  by  weak  convergence  in  L2(0,1),  u(TQ,x),  ux(T0,x),  uxx(Tg,x),  UxxxfTg.x', 

ufc(T0,x),  utx(T0,x) ,  utxx(Tg»x>  are  a11  in  L2 (0,1).  But  then,  using  u(TQ,x), 

ut(T0.x)  as  new  initial  data,  we  may  extend  u(t,x)  to  some  interval 

[T0,T0  +  e] ,  beyond  Tg,  and  this  is  a  contradiction  since  [0,Tg)  is  assumed 

maximal.  The  function  u(t,x)  will  be  a  solution  of  (2.7),  with  w(t,x)  =  u(t,x), 

and  thus,  as  noted  above,  u  (t,*),  u^^  (t,»),  u  (t,*),  u  (t,*)  are  all  in 

ttt  ttx  txx  xxx 

C([0,T]»  L2 ( 0 , 1 ) ) ,  for  every  T  in  (0,Tn).  The  proof  is  complete. 


3.  Global  Solutions.  Our  objective  in  this  section  is  to  show  that  when  the  body 
force  is  "small"  the  maximal  interval  of  existence  of  solution  to  (1.9),  (1.11), 
(1.18)  is  (-00  ,  00  )  and  solutions  decay  as  t  -*■  ®.  For  that  purpose,  the  dissipative 
character  of  viscosity,  embodied  in  assumptions  (1.3)  on  the  relaxation  function  i , 
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plays  the  crucial  role.  Assumption  (1.3)  will  be  exploited  here  through  its 
consequences  recorded  in  the  following 

Lemma  3.1.  There  exist  positive  constants  6 ,  y  such  that 

st  2  s  t 

(3.1)  /  [/  a(t  -  T)w(T)dx]  dt  <  B  f  w(t)  /  a(t  -  i)w(T)dtdt  , 

■00  —00  ■OO  — oo 

St  2  st 

(3.2)  /  [/  a' (t  -  T)w(T)dT]  dt  <  y  f  w(t)  /  a(t  -  T)w(T)dTdt  , 

—00  — OO  —00  —00 

2 

for  any  s  e  (-*,*)  and  every  w(t)  e  L  (-“>,s). 

The  proof  can  be  read  off,  for  example,  from  Lemma  4.2  of  [7],  recalling  that 
a(t),  a'(t),  a"(t)  are  in  L^O,*),  and  that  (by  assumption  (1.3)) 
a(t)  -  a  exp(-t)  is  a  positive  definite  kernel  on  [0,“)  for  some  a  >  0.  As  a 
matter  of  fact,  we  may  use 


(3.3) 


(3.4) 


8  *  —  {/  |a(t)|dt}2  +  “  {/  |a'(t)|dt}2  , 

“  0  °  0 

00  00 

Y  =  ~  {/  |a'(t)|dt}2  +  ~  {/  |a"(t)|dt}2  . 

0  0 


Another  important  implication  of  a  combination  of  (1.3),  (1.5)  and  (1.6)  is  the 
property : 

Lemma  3.2.  Let  k(t)  be  the  resolvent  kernel  of  the  operator 

t 

(3.5)  '  (O)oi(t)  +  /  a*(t  -  t)<*’ (0)u(T)dT  s 

—00 

that  is,  k  is  the  unique  solution  of  the  linear  Volterra  equation 

t 

(3.6)  <p  *  (0  )k(  t)  +  /  a'(t  -  T)ip' (0)k(T)dt  =  -^'(O)a'(t)  . 

0 

Then  k(t)  f 
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The  proof  of  Lemma  3.2  follows  by  a  standard  argument:  Since 
a'(t)  e  L^  ( 0  ,*) ,  the  Paley-Wiener  theorem  states  that  k(t)  e  L^O,®)  if  ^nd  only 
if 

def  *  ' 

(3.7)  P(z)  =  y-'(0)  +  *'(0)a'(z)  =  x*(0)  +  ip’  (0 )  za  ( z ) 

does  not  vanish  on  the  half  plane  Rez  >0.  (In  (3.7)  z=  £+i£  and  A  denotes 
the  Laplace  transform) . 

A  simple  calculation  yields 

A  A 

(3.8)  ReP(z)  =  x'(0)  +  ^•(O)CRea(z)  -  (0  )CIma(  z)  , 

(3.9)  IiWP(z)  =  <ji'(0)?Rea(z)  +  i|>' (0  )£Ima(  z)  . 

On  account  of  (1.3),  (1.5)  and  (1.6), 

(3.10)  ReP(C  +  i0)  =  x'<0)  +  ^'(O)Ca(C)  >0,  0  <  £  <  •  . 

A 

As  regards  ImP(z),  since  by  the  strong  positivity  of  a(t),  Rea(i£)  >0,  we  have 

A 

I  trip (0  +  i?)  *  \|>' (0  KRea(  i£ )  is  positive  for  d  >  0  and  negative  for  C  <  0.  On 

the  other  hand,  ImP(£  +  iO  )  =  0,  0  <  5  <  “>.  Furthermore,  since  a'(t)  e  L^O,*), 

we  deduce  by  the  Riemann-Lebesgue  lemma  that  lim  ImP(z) 

1*1-** 

=  ip'(0)  lim  Ima'(z)  =  0,  uniformly  on  Rez  >  0.  Rut  ImP(z)  is  harmonic  on 

|z|+“ 

Rez  >0  so  that,  by  the  maximum  principle,  we  conclude  that  ImP(z)  >  0  on 

{z  *  5  +  Hi?  >  0,  5  >  0}  and  ImP(z)  <0  on  {z  =  5  +  i tU  >  0,  c  <  0}.  In 

conjunction  with  (3.10)  this  yields  P(z)  #0  on  Rez  >  0  and  the  proof  of  the 
lemma  is  complete. 

Before  proceeding  to  the  proof  of  Theorem  1.1,  let  us  show  that  the  boundary 
conditions  (1.18)  are  equivalent  to 

(3.11)  u  (t,0)  =  u  (t,1)  =  0,  -®  <  t  <  *  . 

x  x 

We  multiply  (1.8)  by  if/(e(t,x))t,  integrate  over  (-®,s),  -00  <  s  <  00 ,  and  use  the 
positivity  of  a(t)  to  get 

s 

(3.12)  /  X<e(t,x)  )i>(e(t,x) )  dt  <0,  -®  <  s  <  ”,  x  =  0,1  , 


or 
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(3.13) 


4'(e(  s,x) )  <  0, 


<  s  <  ”,  x  =  0,1  , 


where 


(3.14) 


,  ,  C 

Vie)  ®  /  x(-n)4»*  (n)dn 


On  account  of  (1.5),  (1.6),  'i'(e)  >  0  on  (-6,6)\{0},  <$  positive  small.  Thus 
(3.13)  yields  e(s,x)  =0,  -“  <  s  <  ”,  x  =  0,1,  and  (3.11)  has  been  established. 
Proof  of  Theorem  1.1.  By  virtue  of  (1.4),  (1.5),  (1.6),  there  are  positive  6  and 
<  such  that 


(3.15) 


vJ’(e)  >  tc,  ^‘(e)  >  k,  x'(®)  >  K>  l®l  <  6  . 


We  modify  ^(e)  outside  the  interval  t-6,6]  so  that  (2.2)  be  satisfied  and  we 
let  u(t,x)  be  the  solution  to  (1.9),  (1.11),  (3.11)  on  a  maximal  time  interval 
(-*,T0). 

For  T  e  (  ""“/Tq  ) ,  we  set 

1  2  2  2  2  2 

(3  16)  U(T)  =  sup  /  (t,x)+u  (t,x)+u  (t,x)+u  (t,x)+u  (t,x) 

(-”,T]  0  t  x  tt  tx 

+  UL(t'x)+uttt(t'x)+uttx(t'x,+utxx(t'x)+Uxxx(t'x,}dx 


t  (  2222222  2  2  2 

+  j  /  {u  +u  +u  +u  +u  +u  +u_  +U.  +u  +u  }dxdt  . 

M  Q  t  X  tt  tx  XX  ttt  ttx  txx  XXX 

Our  strategy  is  to  show  that  there  are  positive  constants  v,  R,  v  <  <5,  such  that. 


(3.17)  |ux(t,x)|2  +  |utx(t,x)|2  +  |uxx(t,x)|2  <  v2,  -«  <  t  <  T,  0  <  x  <  1, 


(3.18) 


U(T)  <  KG 


where  G  is  defined  by  (1.23).  Once  this  claim  has  been  established,  we  may 
complete  the  proof  of  the  theorem  by  the  following  line  of  argument  similar  to  that 
previously  used  in  [6]:  First  we  note  that,  by  virtue  of  our  assumptions  on 
v(t,x),  (3.17)  is  automatically  satisfied,  as  a  strict  inequality,  when  t  is 


sufficiently  small.  Next  we  observe  that,  in  view  of  the  Poincare  inequality 
|u  (t,x)|2  +  I u  (t,x) | 2  +  I u  ( t ,  x ) I  2 

a  L  X  X  X 

(3.19) 

<  /  tluxx(t,y)|2  +  futxjc<t,y)  | 2  +  I2)dy  , 

2  2  2 

and  when  G  <  y  with  y  <  v  /K,  (3.18)  implies  (3.17)  (as  a  strict  inequality). 

2  2 

Thus,  for  G  <  y  <  v  /K,  (3.17)  and  (3.18)  will  hold  for  every  T  on  the  maximal 

interval  of  existence  in  which  case  Theorem  2.1  (in  particular  (2.4))  implies 

Tq  =  ”.  From  (3.16),  (3.18)  we  have  u(t,»),  ut<t,*),  ux<t,»),  u  (t,»), 

u.  (t,*),  u _ (t,*),  u  (t,»),  u  (t,»),  u  (t,*),  u  (t,»)  in  L2  ( (-»,») ; 

tx  XX  ttt  ttx  txx  XXX 

L2(0,1))  OL  ((-«,“);  L2(0,1)).  Now  u(t,«),  Ut(t,»),  ux(t,»),  ufct(t,*),  utx(t,»), 
u  (t,*),  u  ( t , • ) ,  u  (t,»),  u  (t,*)  in  L2((-“,”);  L2(0,1))  implies 

XX  ULU  buX  lXX 

2 

(3.20)  u(t,«),U  (t,*),u  (t,»),u  (t,»),U  (t,»),u  (t,»)  0 ,t  ”, 

t  X  tt  tx  XX 

which,  in  conjunction  with  u  (t,«),  u  (t,»),  u  (t,»),  u  ^  (t,*),  u^  (t,*), 

x  tx  xx  ttx  txx 

uxxx(t,,)  in  (-",<»);  L2( 0,1)),  yields  (1.25). 

It  thus  remains  to  verify  (3.18)  under  the  assumption  (3.17).  We  fix  s  in 
( -”,T] •  The  first  estimate  is  obtained  by  multiplying  (1.10)  by  ’^ux^tx' 
integrating  over  (-”,s]  x  [0,1]  and  integrating  by  parts.  The  reader  should  be 
cautioned  that  in  these  and  the  many  integrations  by  parts  which  follow,  there  are 
several  possible  ways  to  carry  out  such  integrations.  The  ones  selected  in  this 
section  are  chosen  for  the  purpose  of  using  the  same  estimates  when  considering  the 
boundary  conditions  (4.1),  (4.2)  below  (see  Theorem  4.1),  The  result  of  this 
calculation  is 
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l 


(3.21)  j/  ♦'  («x<s.x)  )utx(s,x)dx  +  —  /  x*  («x(  s,x) )  ijj '  (u^  (  s,x)  )u^(  s  ,x)  dx 

si  t 

+  /  /  <|»(u  )  /  a(t-T)ip(u)  dtdxdt 

—  o  x  tx  i.  x  Tx 

1  s  1  si 

mjJ  f  *"<ux>utxdxdt  -  j  /  /  {X'<ux)ir(u  ) 

-«>  0  -»  0 

,  1 

-  x^u^*  (ux)}utxuxxdxdt  -  /  g(s,x)\Hux(s,x))xdx 
s  1 

+  /  /  g^i|»(u  )  dxdt  . 

'  t  xx 

••  o 

To  motivate  our  next  estimate,  we  differentiate  (1.10)  with  respect  to  t  and 
then  integrate  formally  by  parts  to  get 

(3.22)  uttt  =  X<«x>tx  +  /  a{t  -  T>*<ttx>TTxdT  +  gt  * 

— 00 

We  would  like  to  multiply  (3.22)  by  <]>(u  )  and  then  integrate  over 

x  ttx 

(-®°,s]  x  [0,1]  in  order  to  arrive  at  an  estimate  analogous  to  (3.21). 
Unfortunately,  this  operation  is  not  legitimate  since  xj»( ux)  does  not 
necessarily  exist  as  a  function.  Consequently,  same  as  with  the  derivation  of 
(2.12)  in  Section  2,  we  shall  have  to  work  first  with  a  discrete  analog  of  (3.22) 
and  then  pass  to  the  limit.  To  this  end  we  apply  to  (1.10)  the  forward  difference 
operator  A,  of  step  n  >  0,  thus  arriving  at 

t 

(3.23)  Au  =  Ax(u  )  +  /  a(t  -  t)A«J»(u  )  dT  +  Ag  . 

tt  X  X  '  X  TX 

—00 

We  now  multiply  (3.23)  by  ^^ux\x'  we  integrate  over  (-®,s]  x  [0,1],  we  perform 
a  number  of  integrations  by  parts,  we  divide  through  by  n“,  and  we  pass  to  the 
limit  as  n  *  0.  The  outcome  of  this  tedious  but  straightforward  calculation  is 
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1  1  2  1  1 
-  /  0'(U  (s,x))u  (s,x)dx  +  -  /  x'<u  (s,x))u~  (s,x)dx 

tuX  ^  q  X  X  txx 


(3.24)  2  q  x 


s  1 


+  lim  —  /  /  At|>(u  )  /  a(t  -  T)Ai);(u  )  dTdxdt 

*  a  X  tX  _  X  TX 


1 

,im  — 

n+0  rf  -*  0 
s  1 


=  2  /  /  <{'"(ux,utxuttxdxdt  '  ^  >T<Ujs,x))ufjs,x)U1_(.Js,x)dx 


-»  0 


tX  ttx 


Si  3  si 

+  /  /  *-(ux)u  U  dxdt  +  f  /  /  (x"(u  )*'(u  ) 

-00  o  X  tx  ttx  ^  -oo  0 


-  x'(u  W'tu  )}u  u.  dxdt 
x  x  tx  txx 


s  1 

+  {  /  {x"(u  )i)i'(u  )  -  x'(u  )i|)"(u  )}u  u  u  dxdt 

g  XX  X  X  XX  ttx  txx 


s  1 

/ 

0 


+  J  /  (x"'  (“x)f  (ttx)-X"(«x)*"(ttx)-X'  (“x)'*'’"  (ux)}utxUxxUtxxdxdt 


s  1 

/ 

0 


S  1 

/  /  X-(-.x)0”<ux)utito^ottxO«at  -  /  /  X " ( «x ) 0" ’  ( ox ) utxuxIdxdt 

— 00  o 


1 

-  /  X"(u  (s,x)  )|J|’  (u  (s,x))u  (s,x)u  (s,x)u  (s,x)dx 
Q  x  x  tx  XX  txx 


1  si 

-  /  g1t(s,x)*(ux(s,x»txdx  ♦  /  J  g1tt^(ux)txdxdt 
U  — 06  o 


s  1 

+  i  l  g2tx*(ux,ttdxdt  * 
-*C0  0 


The  reader  should  note  that  because  iMu  )  does  not  exist  as  a  function,  one 

x  ttx 

cannot,  after  dividing  (3.23)  by  r\ ,  pass  to  the  limit  as  n  -*•  n  under  the 
integral  on  the  right  hand  side  of  (3.23).  However,  the  limit  of  every  other  term 
in  (3.23)  exists  as  n  +  0,  and  therefore, 
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t 

lim  /  a(t  -  r)Ai|/(u  )  (T,x)dr 
n+O  -  x  Tx 

exists  for  t  e  [0,T],  0  <  x  <  1.  The  same  comment  (arrived  at  by  the  same 
reasoning)  applies  to  the  limit  as  n  -*•  0  of  the  multiple  (quadratic  form)  integral 
on  the  left  hand  side  of  (3.24).  It  is  important  for  the  subsequent  estimates  (in 
particular  (3.27),  (3.28))  to  know  that  this  limit  exists  and  is  finite  (in  fact 
positive).  This  is  important  for  the  concluding  part  of  the  proof  of  Theorem  1.1 
(see  argument  preceding  (3.37)  below). 

To  get  our  next  estimate  we  multiply  by  utxx  the  identity 

t  t 

(3.25)  a(0)Ai|i(u  )  =  -  /  a'(t  -  t)AiJ;(u  )  dT  +  /  a(t  -  t)A\J/(u  )  dt 

XX  1  XX  1  X  TX 


and  we  integrate  over  (-»,s]  x  [0,1].  We  majorize  the  right-hand  side  of  the 
resulting  equation  by  first  applying  Schwarz's  inequality  and  then  using  (3.1)  and 
(3.2).  The  result  is 

(3.26)  a(0)  /  l  UtxxA*(Ux,xdxat 

— flo  o 


</  I 

-at  0 


u^  dxdt}  ^ 
txx 


w  I  I 

-at  0 


A\J»(  u  ) 

XX 


a(t-T  )A<J>(ux)^dTdxdt} 


s  1 

<f  / 

-at  0 


dxdt} 

txx 


v2 


s  1 

(6  /  / 

-at  0 


A*{ux,tx  1 


a(t-T)AiJ)(u  )  dtdxdt} 

X  TX 


v2 


Dividing  through  by  n  (the  step  of  the  forward  difference  operator  A),  letting 
rj  +  0,  and  using  (3.15)  and  the  Cauchy-Schwarz  inequality,  we  end  up  with  the 
estimate 
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(3.27)  J<a(0)/  l  UtxxdXdt'^m/  l  ♦‘Vtx  J  a<t-T)*(u  )  dxdxdt 

— no  Q  —oo  0  — » 


S  1 


<a(0 


7  lim  "^T  /  /  A,J,(ux)tx  ^  a< fc-T  ) Aip(ux) Txdxdx:c3t 


n+o  H  -*  0 


s  1 

<  -  a(0)  /  /  i|/" (u  )u  u  u  dxdt  . 

'  i  X  tx  XX  txx 

-00  0 


Next  we  integrate  over  (-“>,sj  x  [0,1]  the  square  of  (3.23),  we  use  (3.1), 


then  we  divide  through  by  n  and  we  let  n  *  0  to  get 


s  1 


s  1 


/  J  <ttdxdt-4/  /  X’<u)2u2  dxdt 

-“  0  -“  0 


(3.28) 


s  1  t 

-  40  lim—/  /  A\|i(uv)+.v  /  a(t-i)AiJi(ux)TxdTdxdt 


2  j  j  “ r '“x ' tx 

n+o  n  -»  o  x  tx 


Si  55?  S  1 

<  4  /  /  x"(u  )  u  u  dxdt  +4  J  /  g  dxdt  . 

-c  o  x  tx  xx  i.  0 


To  the  above  estimate  we  append 


s  1  _  si 

(3.29)  /  J  u  dxdt  -  /  /  u 

-00  0  -00  0 


which  can  be  derived  from 


(3.30) 


s  1  _  si 

/  /  Autxdxdt  =  /  / 


-00  0 


by  passing  to  the  limit.  In  turn, 
parts. 


1 


.  . . u.  dxdt  +  / 
ttt  txx  ^ 


( s ,x)u  (s,x)dx  =  0 
tt  txx 


1 

Au  Au  dxdt  -  /  Au  (s,x)Au  (s,x)dx 
tt  xx  i.  t  xx 


(3.30)  can  be  easily  verified  via  integrations  by 
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We  now  differentiate  (1.9)  with  respect  to  t 


(3.31) 


u 


ttt 


*(Vtx  +  !  a’(t  -  T)*(ux>TxdT  +  9t  ' 
*00 


and  we  easily  get  the  estimate 


(3.32) 


2  2 
/  uttt(s,x)dx  -5  /  >f '  (uv(s,x)  )uVvw(  s,x)dx 


txx 


-  5 {/  |a’(t)|dt}  sup  /  i|>'(u  )u2  dx 

*  Y  tYY 

0  <--,3]  0 


1 


<  5  /  ^"(u  (s,x) )u  (s,x)u  (s,x)dx 

Q  X  tx  XX 


+  5 {/  |a'(t)|dt}2  sup  /  iji"(u  )2U2  U2  dx  •*  /  g2(s,x)dx  . 

0  0  *  tx  xx  o  t 


The  final  set  of  estimates  is  derived  by  the  following  procedure:  we 
differentiate  Equation  (1.9)  with  respect  to  x  and  then  add  and  subtract 
appropriate  terms  to  arrive  at 


(3 


<p'(0)u  +  /  a'(t  -  t)<|>'(0)u  dT  «  u  -  [*'  (u  )  -  if*  (0)]u 

•  33)  xxx  ‘  T  xxx  ttx  x 


XXX 


-  *>"(u  )u  -  /  a'(t  -  T)[\J|'(U  )  “♦'(0)]u  dt 
X  XX  1  _  X  xxx 


-  /  a'(t  -  T)ij>"(u  )uz  dT  -  g  *  X(t,x)  . 

_  X  XX  X 


Thus,  if  k(t)  is  the  resolvent  kernel  of  the  operator  (3.5), 


t 

(3.34)  s5'(0)u  x<t,x)  =  X(t,x)  ♦  J  k(t  -  T)X(T,x)dT  . 

— m 

ny  Lemma  3.2,  k(t)  f  so  that  we  have  estimates 


(3.35) 


0  '  O  1  o 

s? '  ( 0  )  /  u  ( s ,  x )  d x  <  2  /  X  ( s ,  x ) dx 

i  xxx  ' 

0  0 

»  1 
+  2 {/  |k(t)|dt}2  sup  /  X2(t,x)dx  , 

0  (-»,s]  0 

?  °  1  2  °°  2  S  1  2 

(3.36)  v?’(0)  /  /  u  dxdt  <  2{l  +  [/  |k(t)|dt]  }  /  /  X  (t,x)dxdt  . 

_oo  o  xxx  0  -»  0 

We  are  now  ready  to  prove  (3.18)  under  assumption  (3.17).  First  we  note  that, 
on  account  of  (3.15),  U(T)  can  be  majorized,  with  the  help  of  the  Poincare 
inequality,  by  the  supremum  over  (-“,T]  of  an  appropriate  linear  combination  of 
the  left-hand  sides  of  the  estimates  (3.21),  (3.24),  (3.27),  (3.28),  (3.29),  (3.32), 
(3.35)  and  (3.36).  On  the  other  hand,  each  term  on  the  right-hand  sides  of  these 
estimates  can  be  majorized,  by  means  of  the  Cauchy-Schwarz  inequality  and  (3.17),  by 
either  cG,  or  o(v)U(T),  or  eu(T)  +  c(e)G  for  any  e  >  0.  We  thus  arrive  at  an 
estimate  of  the  form 

(3.37)  U(T )  4  {0(v)  +  0(e) }u(T)  +  c(e)G 

from  which  one  can  get  (3.18)  by  fixing  v  and  e  sufficiently  small.  The  proof 
of  Theorem  1.1  is  complete. 

4.  Remarks  and  Extensions.  When  the  endpoints  of  the  body  are  pinned,  in  the  place 
of  (1.18)  we  have  boundary  conditions 

(4.1)  u(t,0)  =  u(t,1)  =0,  —  <  t  <  •  . 

Similarly,  when  one  endpoint  (say  x  =  0)  is  pinned  and  the  other  is  free, 

(4.2)  u(t,0)  =  0,  o(t,1)  =0,  <  t  <  •  . 

In  these  cases  no  rigid  motions  are  possible  so  we  don’t  have  to  assume  (1.20)  nor 
should  we  expect  that  (1.21)  will  generally  hold. 

When  the  body  force  satisfies  (1.22)  with  g2(t,x)  =  0,  all  estimates  derived 
in  Section  3  for  the  case  (3.11)  are  also  valid  under  (4.1)  or  (4.2).  We  thus  have 
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Theorem  4.1.  There  is  u  >  0  with  the  property  that  for  every  q(t,x)  on 

(-“,“)  *  [0,1],  which  satisfies  (1.22)  with  g2(t,x)  i  0  and  (1.24),  and  any 
v(t,x)  on  ( ~® , 0 ]  x  [0,1],  with  v(t,*),  vt<t,*),  v^ft,*),  v  (t,»),  v  (t,»), 

<t,*),  v  (t,*),  v  (t,»)  in 
XX  ttt  ttx  txx  XXX  — 

C  ((-<“,  0  ] ;  L2 ( 0 , 1 ) )  n  L2 ( ( -® , 0 ] ;  L2 ( 0 , 1 ) ) ,  which  satisfies  (1.9)  together  with 

(4.1)  ( or  (4.2))  for  t  <  0,  there  exists  a  unique  u(t,x)  on  (-<*>,»)  x  [0,1], 

with  u( t , • ) ,  u  (t, • ) ,  u  (t ,  • ) ,  u  ( t , • ) ,  u  ( t,  •) .  u  (t, • ) .  u  ( t, • ) ,  u  (t. • ) . 
-  t  x  tt  tx  xx  ttt  ttx 

u_(t,«),  u  (t,*)  in  C(  ( ~00,00) ;  L2(0,1))  n  L2  ( ( -00  ,°o ) ;  L2(0,1),  which  satisfies 
(1.9),  (1.11)  and  (4.1)  (or  (4.2)).  Moreover ,  (1.25)  holds . 

As  history  and  body  force  get  smoother,  solutions  become  smoother.  Regularity 
results  can  be  obtained  by  establishing  a  priori  estimates  for  derivatives  of  u  of 
order  4,  5,  etc.  Such  estimates  fall  into  three  categories:  those  derived  by 
differentiating  (1.10)  a  number  of  times  with  respect  to  t  and/or  x  and  then 
multiplying  by  the  appropriate  multiplier  (recall  the  derivation  of  (3.24));  those 
derived  by  expressing  certain  derivatives  in  terms  of  other  derivatives  through  the 
equation  itself  (compare  with  (3.32)  or  (3.35));  those  obtained  through 
interpolation  (such  as  (3.29)).  The  program  is  feasible  because,  since  the  problem 
is  autonomous  (kernel  of  convolution  type),  differentiations  with  respect  to  x 
or  t  essentially  preserve  the  form  of  the  equation  (compare,  for  example,  (3.22) 
with  (1.10)).  Time  derivatives  of  u  satisfy  the  same  homogeneous  boundary 
conditions  as  u,  at  x  *  0,1,  so  that  differentiating  the  equation  with  respect 
to  t  is  generally  a  better  prospect  than  differentiating  with  respect  to  x.  In 
any  event  there  are  so  many  possible  combinations  of  differentiations,  integrations 
by  parts,  etc.,  that  one  may  establish  several  variants  of  regularity  theorems. 

Here  is  a  typical  one: 


Theorem  4.2.  Suppose  the  assumptions  of  Theorem  1.1  hold  and,  in  addition,  $  and 

are  C4  smooth,  v _ (t,»),  v.  .  .  (t,*),  v.  .  (t.»),  v  (t.»)  v  (t.»)  are 

-  -  tttt  tttx  ttxx  txxx  xxxx  - 

in  C(  ( -*,0  ] ;  L2(0,1))  n  l2  ( (  —  ,0] ;  L2(0,1))  and 


(4.3) 


gtt(t,»),gtx(t,*)/gxx(t,»)  in  C((-o»,«)jL  (0,1))^ 

l2((-»,»).  t2(0,1>) 


g(tfX)  =  g1 (t,x)  +  g2(t,x)  with  > '92ttx(t ' *  * 


in  L2(0,1))  . 


Then ,  when  (1.24)  is  satisfied  with  y  sufficiently  small,  the  solution  u(t,x)  of_ 
(1.9),  (1.11),  (1.18)  possesses  u  (t,»),  u  (t,*),  u  (t,*),  u  (t,»), 

t t t »  t v uX  lCXX  uXXX 


Uxxxx(t '  *  ’  —  n  L2((— ,-);  L2(0,1))  and 

(4-4)  uttt(t'*)'uttx(t'*,'utxx<t',)'uxxx(t'*)  iSHr  °'  *  *  •  • 

It  is  noteworthy  that  the  extra  derivatives  (4.3)  of  g  that  are  required  in 

order  to  guarantee  smoothness  of  the  solution  need  not  be  "small".  This  is  due  to 

the  fact  that  all  energy  integrals  are  quadratic  forms  in  the  higher  order 

derivatives  of  u,  with  coefficients  that  are  solely  controlled  by  v  of  (3.17). 

As  we  have  seen  in  Section  3,  v  is  controlled  by  U  which,  in  turn,  is  controlled 

by  G. 

We  close  with  remarks  on  the  multidimensional  situation.  The  configuration  of 
the  body  is  now  a  bounded  set  !)  C  sn  with  smooth  boundary  9ft  and  the 
displacement  is  an  n-dimensional  vector  field  u.  A  typical  problem  is  to  determine 
u(t,x),  <  t  <  *,  x  e  0,  such  that 


(4.5) 


n  9*  t  91 

I  {aT-1  +  /  ■’  (t  -  t)  t—-1  dr}  +  g.  ,  i  «  1 , . . .  ,n  , 

j-1  3  j  -•  3xj  1 

_oo  <  t  <  “,  xeft, 


(4.6) 


(4.7) 


u( t ,x)  »  v( t ,x) ,  ■•  <  t  <  0,  xeft, 


u(t,x)  -  0, 


*•  <  t  <  •,  X  e  9ft  , 


where  4  and  f  are  known,  smooth,  matrix  valued  functions  of  the  matrix  g  *  Vg 
(strain),  g(t,x)  is  an  assigned  body  force  and  £(t,2e)  is  the  given  history. 
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We  assume  4(0)  =  *F(0)  =  0,  set 


(4.8) 


3$.  .(e) 

C.  ..  „(e)  =  -  ,  D.  ...(e) 

i;jkx,  ~  3ej{£  xjkl  ~ 


(4.9)  E....(e)  *  C,  ...(e)  -  a(0  )D.  (e)  , 

i]k£  ~  ljki  ~  i}k£  ~ 

and  impose  the  symmetry  restrictions 


(4.10) 


Cijk£  Ck£ij'  Dijk£  Dk£ij 


Assumptions  (1.4),  (1.5),  (1.6)  will  here  turn  into  coercivity  conditions  for  the 
partial  differential  operators  associated  with  Cijk£^'  Dj.jki/~*  and  Eijk£^~*’ 
Under  boundary  conditions  (4.7),  coercivity  is  equivalent  to  strong  ellipticity 


(4.11) 


.  I  ,  'f  -  's'  * 1 


(4.12) 


,  A, ,  “uM'a’Wi'* > - 1  • 

i,  j  ,k  ,x 


(4.13) 


,  A  .  '‘linWWft  >  'S'  -  If  -  ’  • 

1 , 3  »*  ,* 


Assumptions  (4.10),.  (4.11),  (4.12)  and  (4.13)  can  be  motivated  by  Mechanics. 


However,  in  order  to  carry  through  the  analysis,  we  require  an  additional  condition 
whose  physical  interpretation  is  less  clear.  We  define 


(4.14) 


jkipqr 


(e)  “Id  .(e)E  (e) 

~  7  ilk£  ~  ipqr  ~ 


and  assume  that  F  is  symmetric, 

(4.15)  F...  «  F  , 

jk£pqr  pqrjkf 

and  that  its  value  at  e  =  0  corresponds  to  a  coercive  operator.  We  note  that  in 
the  special  case  4  =  ¥  the  resulting  F  automatically  satifies  the  above 
conditions. 

Under  the  above  assumptions  it  is  possible  to  establish  the  existence  of 
globally  defined  smooth  solutions  to  (4.5),  (4.6),  (4.7)  by  the  procedure  followed 
here  in  the  one-dimensional  case.  The  strateay  is  to  establish  a  priori  eneray 


2 

estimates  for  the  L  (Q)  norms  of  derivatives  of  u  of  sufficiently  high  order 
(depending  upon  n)  that  would  guarantee,  via  Sobolev's  lemma,  pointwise  bounds 
analogous  to  (3.17).  The  calculations,  however,  are  very  long. 

It  is  easy  to  discern  the  essential  ingredients  in  the  proofs  and  it  thus  seems 
feasible  to  develop  an  existence  theory  for  the  history-value  problem  in  a  class  of 
abstract  nonlinear  integrodif f erential  eouations 

d2 

(4.16)  — |  =  A( u( t) )  +  /  a' (t  -  T)B(u(T))dx  +  g(t) 

dt2  -“ 

on  a  Hilbert  space  H,  where  A  and  B  are  nonlinear  operators  defined  on  a  scale 
of  Hilbert  spaces  (abstracting  the  scale  of  Sobolev  spaces  [w  '  (SI)]  )  and 
satisfying  appropriate  symmetry  and  coercivity  conditions.  We  remark  also  that  the 
general  and  physically  interesting  case  in  which  the  stress-strain  relaxation 
function  is  a  n  *  n  matrix  a  (in  (4.5)  a  =  al)  is  considerably  more 
complicated  than  the  situation  considered  here. 
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